Quantum criticality between topological and band insulators in (3 + 1) -dimensions 
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Four-component massive and massless Dirac fermions in the presence of long range Coulomb interaction 
and chemical potential disorder exhibit striking fermionic quantum criticality. For an odd number of flavors of 
Dirac fermions, the sign of the Dirac mass distinguishes the topological and the trivial band insulator phases, 
and the gapless semi-metallic phase corresponds to the quantum critical point that separates the two. Up to 
a critical strength of disorder, the semi-metallic phase remains stable, and the universality class of the direct 
phase transition between two insulating phases is unchanged. Beyond the critical strength of disorder the semi- 
metallic phase undergoes a phase transition into a disorder controlled diffusive metallic phase, and there is no 
longer a direct phase transition between the two types of insulating phases. Our results are also applicable to 
even number of flavors of Dirac fermions, and the band inversion transition in various non-topological narrow 
gap semiconductors. 



The low energy, long wavelength quasi-particle spectrum of 
various narrow gap semiconductors is well approximated by 
noninteracting (3 + 1) -dimensional massive Dirac fermions 
and by adjusting the chemical composition, or by applying 
pressure, the sign of the Dirac mass can be changed at a band- 
inversion transition (BIT). At BIT the system becomes semi- 
metallic (SM) and is described by massless Dirac fermions. 
The SM phase is an interesting example of z = 1 fermionic 
quantum critical point (QCP). In the massless phase, the con- 
duction and valence bands cross at a discrete number of di- 
abolic points inside the Brillouin zone and depending on the 
number of inequivalent diabolic points we obtain multiple fla- 
vors of Dirac fermions. The narrow gap semiconductors such 
as Pbi-ajSn^iTe (four flavors), Bii_j;Sbj; and Hgj^_^Cd2;Te 
(each with single flavor) are well known examples, which for 
special values of x, become massless 111. For odd number 
of Dirac fermion flavors, the QCP describes the phase transi- 
tion between a topological insulator (TI) and an ordinary band 
insulator (BI) ||2|-|5|]. The recent experimental observation of 
TI phase in different narrow gap semiconductors have spurred 
our interest in the (3 + 1) -dimensional Dirac materials US]. 
The low energy spectrum of materials like Bii-j^Sba;, Bi2Te3, 
Bi2Sc3, where TI phase has been observed are all described 
in terms of a single flavor massive Dirac fermion 

The stability of the disorder free SM phase in the presence 
of long range Coulomb interaction and the noninteracting SM 
phase in the presence of random chemical potential disorder 
were respectively addressed in Refs. Jst] and |@]. Recently 
the noninteracting, disordered (3+l)-dimensional TI has been 
considered in Refs. lflo[[Tlll . For the noninteracting disordered 
problem a symmetry based ten-fold classification of TI and su- 
perconductors has been described in Ref. lll2ll . and the stabil- 
ity of the two dimensional surface states has been discussed on 
the basis of this symmetry classification. The noninteracting 
TI in the presence of generic time reversal symmetric disorder 
belongs to the symplectic All class. In general such classifi- 
cation does not hold in the presence of interaction. The sta- 
bility of the two dimensional disordered and interacting sur- 



face states of a TI has been addressed in Ref. |13ll. However, 
the combined effects of interaction and disorder on the bulk 
fermions has not been considered before. Motivated by this 
and possible future experiments in which novel QCPs can be 
explored, we analyze the problem of both massive and mass- 
less Dirac fermions in the presence of Coulomb interaction 
and random chemical potential type disorder using a pertur- 
bative renormalization group (RG) analysis. Remarkably, the 
vanishing density of states at the Dirac points renders such a 
calculation reliable in comparison to the corresponding non- 
relativistic problem. 

For orientation we first consider the disorder free nonin- 
teracting Dirac fermion action. For simplicity we consider 
only one species of Dirac fermion. Assuming inversion (par- 
ity) and time reversal symmetry and using a spinor basis 
ijj^ = (c+,-t', c_.-f-, c_4), where c±,s respectively cor- 
respond to the annihilation operators for parity even and odd 
states, with spin projection s, we can write the following Eu- 
clidean action 



5*0 



d xip 



70 (So - + iv-fjdj + m - Bd^ 



i'- (1) 



The latin index j is a spatial index and A ~ 1/a is the ul- 
traviolet cutoff, where a is the lattice spacing. The parameter 
V is the Fermi velocity and m is the Dirac mass. The anti- 
commuting Euclidean 7 matrices satisfy {7^,7,^} = 2(5^ i^, 
and i/i = i/'^7o. The four-component spinor structure arises 
from the two-sublattice crystallographic structure and the two 
spin components. We have also incorporated two higher gra- 
dient terms involving A and B. The above action is invariant 
with respect to the parity {V) and time reversal (70 transfor- 
mations: Vil}'P~^ = 7o'0, TipT^^ = — 7i73'0- For A = 0, 
the action is also invariant under charge conjugation (particle- 
hole) transformation (C): CtpC^^ = —72V'- The particle-hole 
symmetry breaking term does not affect the topological prop- 
erties, and can be off-set by adjusting of the chemical poten- 
tial, and henceforth we will set ^4 = 0. 

The fermion mass m and the higher gradient term Propor- 
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tional to B break the f/(l) chiral symmetry of the massless 
Dirac fermions defined hy e'^''/^'>'^^il;, ^ ^e'^" 
where 75 ~ 7o7i7273- The TI and BI phases are respectively 
defined by the conditions mB < and mB > and are sep- 
arated by a finite chiral angle 69 ~ tt. This is reflected in the 
quantized magneto-electric coefficients tt and 0, respectively 
for TI and BI. At the critical point m = 0, the dynamic expo- 
nent z = 1, and in the RG sense the higher derivative terms 
can be ignored. 

The Dirac structure of the Hamiltonian allows various types 
of disorder The constraint of time reversal invariance allows 
the following six bilinears: ipjotp, 'iptp, V'TolsV' ™d 'ipjolj'4' 
{ j = 1,2, 3). The bilinears ■ipjoip, ipip respectively correspond 
to random chemical potential and random mass scattering. 
The physical description of other four bilinears depends on the 
crystallographic details. We shall concentrate on the random 
chemical potential as the dominant elastic scattering process, 
and add Sd ~ J d'^xV{x)'tljjo'ip to the action 5*0. The random 
potential V{x) is a Gaussian white noise distribution specified 
by the disorder average ((l/(x)l/(x ))) — Ay(5'^(x— x ). The 
detailed analysis for generic time reversal symmetric disorder 
is provided in the supplementary material il4ll . 

Since typically v/c ^ 10~^ — 10"^ (c is the velocity of 
light), the Coulomb interaction is instantaneous. Its strength 
is characterized by the the dimensionless parameter a = 
e^/(47reu) ~ 2.2/e — 22/e, where e is the static dielectric 
constant of the material. We perform disorder average using 
the replica method, which we use merely as a book-keeping 
device for perturbative RG calculations. The replicated Eu- 
clidean action after disorder averaging of the partition func- 
tion is 



■^a {7o(<9o + ig'Pa) + W7j9i + m - ijja 



X (?/S670'/'6)(,,4)(2) 



where g = y/inva, and a, b are replica indices. We 
have used the abbreviated notation (?/'a7oV'a)(^ ^ ^ = 

^a{x,xo)'jo^a{x,xo)- We have introduced an auxiliary 
scalar potential (pa to decouple the four-fermion Coulomb in- 
teraction term. Unlike two dimensions, the Coulomb interac- 
tion manifests itself as (djifa)'^, which is analytic in momen- 
tum and g can certainly receive loop corrections. The action 
S preserves all three discrete symmetries V, T and C. For 
RG calculations we replace the couplings by the correspond- 
ing dimensionless couplings m ^ m/ (vA), B — > BA/v and 
Av AvA/{2tt'^v'^). The details of the RG calculation are 
provided in the supplementary material lfl4ll . 

The density of states for massless and massive problems 
are, respectively, p{E) oc E"^ and p{E) oc \E\\/ E'^ — m?. 
Since the the density of states vanishes at zero energy, the 
scattering rate t~-^{E) calculated from lowest order Born ap- 
proximation also vanishes at zero energy. At the tree level. 



B and Ay are irrelevant couplings, and a and m are respec- 
tively marginal and relevant couplings. Therefore we antic- 
ipate that the universality class of the QCP between TI and 
BI will be unchanged up to a critical strength of the disorder. 
This should be contrasted to the two dimensional problem, 
where the chemical potential disorder is a marginally relevant 
perturbation and invalidates the lowest order Born approxima- 
tion result I15II . 
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FIG. 1: (a)The RG flow and (b) the phase diagram in the a — Ay 
plane for massless Dirac fermions with m = B = 0. The CDM 
phase represents a disorder induced massless, compressible diffusive 
metallic state with finite density of states and scattering rate at zero 
energy. 

Consider first the massless Dirac fermions with 5 = 0. 
From a momentum shell renormalization group (RG) calcula- 
tion to 0{c?), 0{A\), and 0{Aya), we find 



dv 
H 
da 
H 

dAv 
dl 



2a 
3^ 



, ^ 4a 
— = alAv- — 



Av{-l + 2Ay-^Y 



(3) 
(4) 
(5) 
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By keeping v fixed we obtain a scale dependent dynamic 
exponent z{l) = 1 + Av{l) — 2a{l)/3TT. There are two 
fixed points: (i) attractive, noninteracting, clean fixed point: 
Ay = a = 0, z = 1; and (ii)noninteracting finite disorder 
critical point: Ay = 1/2, a = , z = 3/2. For the nonin- 
teracting disordered problem, the fixed point (ii) controls the 
transition between the two phases where disorder is respec- 
tively irrelevant (SM) and relevant. In the phase where disor- 
der is relevant, both the zero energy density of states, and the 
zero energy scattering rates are finite. Therefore this disorder 
induced phase will be termed a compressible diffusive metal 
(CDM). The RG flow and the associated phase diagram are 
respectively shown in Fig. |l(a)| and Fig. |l(b)| 

By linearizing the flow equations in the vicinity of the fixed 
point (ii), we findQ:(/) w aoe'^^, and Ay - 1/2 - 8a/(37r) w 
(Ayo — 1/2 — 8ao/(37r))e'. Therefore at the disorder con- 
trolled critical point a is a relevant perturbation, and a shifts 
the SM-CDM phase boundary to a larger value of Ay = 
1/2 + 8Q;/(37r). Along this phase boundary the correlation 



length diverges with an exponent v = \, but the dynamic 
exponent varies continuously as z = 3/2 + 2a/7r. Notice 
that CDM phase is a strongly interacting state of matter, and 
this was not addressed in Ref. Igl-fTll In the SM phase, the 
Coulomb interaction initially grows before curling back to- 
wards zero in a logarithmic manner, and the initial growth of 
a is controlled by the bare strength of Ay. This unusual flow 
will be reflected as a non-monotonic temperature dependence 
of the inelastic scattering rate. The critical behavior at the SM- 
CDM phase boundary should be contrasted with its (2 + 1)- 
dimensional counterpart. In (2 + l)-dimensions, there is no 
perturbative loop correction to g, and the phase boundary is a 
line of critical points with z ~ 1^^- 

Now we consider the role of m and B. Compared to the 
SM phase, we expect TI and BI to be stable up to a larger 
disorder, Ay(m) > (1/2 + 8a/(37r)). Due to the irrelevant 
nature of B, we expect it to cause non-universal shift of the 
phase boundaries, leaving the critical properties unchanged. 
For finite m and B, the RG equations are given by 



dv 
ll ' 

dm 
HT 

dB 

IT 

da 
'dl ' 

dA^ 



z - 1 



2a 



A. 



37r^l + (m + B)2 I + [m + B) 

a Ay 
37r^l + (m + B)2 1 + (m + Bf 



B 



a 



7r^l + (m + S)2 1 + (m + Bf 



-B 



37r^l + (m + Bf 
Ay 



37r^l + {m + Bf 
2a 2al + |(m2 + S2) + mB 



l + {m + BY i-K^l + {m + Bf Stt [l + (m + B)2]5 

2Ay 4a 'ial + Um'^ + B^) +mB 



-1 + 



I + {m + Bf 37r^l + {m + Bf Stt [i + (m + 5)2]l 



(6) 
(7) 
(8) 
(9) 
(10) 



In Fig. 2(a)| we show the phase diagram for B = a = Q and 
in Fig. |2(b) we show the phase diagram for a = and a bare 
value i?o = 0.5. When disorder is irrelevant, there is a di- 
rect phase transition between TI and BI phases along the line 
ah. In this region, Ay(Z) ^ Aye^', and B{1) ~ Boe^\ 
and a{l) ^ ao{l + 4ao//(37r))^^, and only relevant vari- 
able is fermion mass m. In this region for ao = 0, we find 
{m{l) - Ay(/)B(/)/3) w eJ{mo - AyoBo/3). Therefore 
(to — Ay_B/3) behaves as the effective mass, and for finite 
B, the TI-BI phase boundary shifts to m = BAy /3. This 
can be seen by comparing the segment ab in Fig. |2(a)| and 
Fig. |2(b)| The Coulomb interaction causes additional shift 
to TO = i?Ay/3 + Ba/{2TT). Therefore in the weak disor- 
der regime, it is possible to induce a transition between two 
insulating phases by tuning the strength of the disorder As 
Ay and a are respectively irrelevant and marginally irrele- 



vant couplings, z asymptotically approaches unity, and the 
universality class is described by the massless Dirac fermions. 
There are logarithmic corrections to the scaling properties due 
to marginally irrelevant nature of a and it is captured through 
the scale dependent z, and also by a factor e'(ao/a)^/^ for 
the scaling dimension of to. The point is a multi-critical 
point at which the massless SM phase undergoes a transition 
either into CDM or one of the two insulating phases. When 
disorder exceeds the critical strength corresponding to b, there 
is no longer a direct transition between two insulating phases. 
Along the TI-CDM and BI-CDM phase-boundaries respec- 
tively denoted by bd and be, z is non-universal, but the mean 
free path still has the exponent unity. The dashed line be sepa- 
rates the disorder controlled CDM phase into two regions with 
negative and positive effective masses, which do not have any 
physical distinction. For a special case of chiral symmetric 
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FIG. 2: The phase diagrams in the m — Av plane for (a)_Bo = Qo = 
0, (b)_Bo = 0.5, a = 0. The direct transition between TI and BI 
along ah is governed by massless Dirac fermions. When disorder ex- 
ceeds a critical strength, the insulating phases are separated by CDM, 
and transitions along fed, and fee have non-universal critical proper- 
ties. The dashed line fee describes a cross-over between two regions 
of CDM with negative and positive effective masses. 



disorder and _B = 0, the chiral symmetric diffusive metallic 
phase along he becomes distinct from the rest of the CDM 
phase. In that case, due to the presence of additional diffu- 
sive modes, the weak anti-localization correction for a chiral 
symmetric CDM becomes two times larger than that of the 
chiral symmetry breaking CDM ifioll . However for B 7^ 0, or 
generic disorder such distinction is lost and he corresponds to 
a cross-over line. The Coulomb interaction shifts the point h 
to higher strength of disorder, and leads to additional non- 
universal shifts of the phase boundaries, and non-universal 
change of z. 

In the weak disorder regime, the transition between TI and 
BI will be accompanied by interesting critical properties of 
massless Dirac fermions |14]. Since (3 + l)-dimensions is 
marginal, some care is necessary to disentangle slow loga- 
rithmic corrections in many physical quantities. For example, 
the specific heat Cy and compressibility k, instead of being 
proportional to /v^ will be proportional to (r^/u^){l + 
4Qo/(37r) log(uA/T)}~^/^. Similar logarithmic corrections 
in (2 + l)-dimensions have been discussed in Ref. iflill . In 
the high temperature limit, the diamagnetic susceptibility in- 
stead of being a constant has logarithmic enhancement x ~ 



— e^z;/(247r^) log(uA/T). A similar logarithmic correction 
proportional to log(i?) appears in the strong field limit. How- 
ever a finite particle-hole symmetry breaking term A and /i 
will lead to conventional diffusive Fermi liquid behavior in 
the low temperature limit specified by /i/T ^ 1. Therefore 
the critical behavior will be limited to T S> /i. However by 
a careful adjustment of /i, the critical properties can be found 
even in the low T limit. The critical behavior will be found 
even for the massive fermions provided that T > m. In the 
critical regime the inelastic scattering rate ~ a^T is larger 
than elastic scattering rate, and dominates the transport in the 
collision dominated regime tu ^ a^T. A quantum Boltzmann 
equation leads to the conductivity 



30.46T 

alog(l/a) 







1 " 









26.67 



n -1 



log(l/a) 



(11) 



The disorder induced initial growth of a will lead to non- 
monotonic temperature dependence of the inelastic scattering 
rate and the conductivity. 

Our results are obtained for inversion symmetric systems, 
and do not apply for Hgj_^Cd2:Te due to broken inversion 
symmetry and the presence of additional gapless quadratic 
band at the F point. In the presence of inversion symmetry 
breaking Dirac mass such as 7075 and 7o7j, even for a clean 
system one can find a metallic phase, if the inversion symme- 
try breaking mass exceeds m. In the supplementary material 
we have considered the effects of inversion symmetry break- 
ing disorder, and the generic phase diagrams shown in Fig. 2 
remain qualitatively unchanged 11411 . In the strong disorder 
limit we have not accounted for the localization corrections 
for low energy diffusive modes, and such corrections can play 
important role in determining the more accurate scaling be- 
havior in the strong disorder limit. The localization correc- 
tions are expected to drive a further phase transition from the 
CDM phase to disorder controlled insulating phase. The nu- 
merical work in Ref. [ijj and Ref. [l7| for noninteracting prob- 
lem in the strong disorder limit have showed the existence of 
a disorder induced topological Anderson insulator phase. Our 
work suggests that, akin to the conventional metal-insulator 
transition problem lIlSll . the interaction effects become strong 
in the diffusive metallic phase. The effects of strong inter- 
action on metal-insulator transition and topological Anderson 
insulator will be addressed in a future publication. 
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Supplementary Material for EPAPS 

I. RENORMALIZATION GROUP CALCULATION FOR GENERIC TIME REVERSAL SYMMETRIC DISORDER 

In this section we provide the details of the RG calculations for Coulomb interacting Dirac fermions in the presence of generic 
time reversal symmetric disorder. We begin with the disorder part of the action given by 



Sd ~ / d'^x^p 



K)(x)7o + M(x)l4x4 + Vb5(x)7o75 + V"oj(x)7o7i 



(12) 



The disorder bilenears ipjoip and ipip always correspond to random chemical potential and random mass scattering. But, the 
physical significance of the other four disorder bilinears depends on the crystallographic details of the material. When the 
point group symmetry is such that, the Dirac 7 matrices correspond to the Dirac representation, 'ijj'yo'yj^ corresponds to random 
spin orbit scattering. For this reason we will denote the coupling constant for Vbj (x) 7/170 7^^/' as Aso- In the following we 
will demonstrate that all the bilinears are mutually coupled, and terms like tp^o'jjtp will be generated due to the interplay of 
random chemical potential and random mass disorder even if its bare strength was zero. The symmetry properties of the disorder 
potentials are summarized in the Table H] We assume Gaussian white noise distributions for all forms of disorder, and the 
appropriate distributions functions are also summarized in the TableU] 



TABLE I: The disorder distributions and their symmetries. The double angular brackets denote disorder average with respect to Gaussian 
white noise with zero mean. 



Bilinear 


V 


T 


C 




Disorder 




+ 


+ 




/ 


({l^,(x)l^o(x')» = Av5»(x-x') 






+ 


+ 


/ 


((K)5(x)1/05(x'))) = Ao55'(x-x') 




+ 


+ 


+ 


X 


((Af(x)Af(x')» = Am<5'(x-x') 






+ 




X 


((K)«(x)V^o,(x'))) = Aso5,,5'(x - x') 



After disorder average using replica method, we obtain the following Euclidean action S, 



S 



ipa {7o(9o + igifa) + v-fjdj +rn- Bd'^}ipa + ^{djfaf 



Av 
2 



d^xdxodxQ ('0a7O^/'a)(^ ,,^) {'ipb'loA)(^^^ 



/a f I - - 

/ d^xdxodx'o (ViaV'a)(,_,^) i^b^b)^^^^ 



d^xdxodXQ (-0a7o75V'a)(^_^^) (V'b7o75'0b) (^^^^) 
d^xdxodXQ (V'a707jV'a)(^^^^) i'^blOjj'ipb)^^^^) (13) 

Now we perform the one loop RG calculation in the momentum shell scheme, where we eliminate the fast degrees of freedom 
belonging to the shell Ae~' < |k| < A, —00 < fcg < oo- All the relevant Feynman diagrams which do not vanish in the replica 
limit n — > are shown in Fig. [3] 

The clean, noninteracting fermion propagator (^□(/co, k) and the bare scalar potential propagator Do(fco! k) are given by 



Go(fco,k) 



-i(7ofco + v^'jkj) + [m + Bk'^) 



kl + w2k2 + {m + Bk 



2\2 



Do(fco,k) = ^ 



(14) 
(15) 
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Disorder induced fermion self-energy : We first consider the disorder induced fermion self-energy diagrams of Fig. |3la). 
Due to the rotational invariance, the disorder induced self-energy diagrams are independent of external momentum and depend 
only on the external frequency. The expression for the disorder induced self-energy is given by 



I]^(fco, k) = ^ A„ y ^r„Go(fco, q)ra 
(Ay + Am + Ac 



3Aso)A/ . (Ay 



3A5o)A/ 



(27r2-y2)(l + {m + BA2)2t;-2A-2) 



(27r2w2)(i + (m + SA2)2i;-2A-2) 



^^^(m + i?A2), 



(16) 



(17) 



where / (Pq/{2TiY denotes the momentum shell integration with Ae ' < g < A. The two parts of I]^(fco,k) respectively 
cause field and mass renormalizations due to disorder. 

Coulomb interaction induced fermion self-energy : Now consider the Coulomb interaction induced exchange self-energy 
graph shown in Fig.[3lb). Due to the instantaneous approximation for the Coulomb interaction this diagram is independent of 
the external frequency. The expression for the exchange self-energy is given by 



S-(fco,k) = V r I 7^7oGo(go,q)7oi?o(fco-9o,k-q) 
J-oo (27r) J {2ttY 



'9 



Up to the quadratic order in k, we obtain 

E^-(fco,k)«- 



27rv/l + (m + SA2)2x;-2A- 



(27r)3 

2^ J (27r)3 [ql + w2q2 + (m + Bq^)] (k - q)2 



4 2 

i-v-fjkj + 2(m + BA^) + -(m + BK^)k-'^k'^ 
3 3 



(18) 
(19) 



(20) 



The disorder and instantaneous Coulomb interaction destroy the Lorentz invariance of the noninteracting action, and this 
is manifested in the inequality of the coefficients of 7ofco and 7jfcj. Together these terms give rise to a multiplicative field 
renormalization constant Z^,, and a scale dependent dynamic scaling exponent z(l). Since the field renormalization at one loop 
order arises only from the disorder contribution. Ward identity ensures that the renormalization of the Coulomb interaction vertex 
comes solely from the disorder contribution, shown in Fig.[3ji). 
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Scalar potential self-energy: We next consider the scalar potential self-energy diagram shown in Fig. [3jc). The scalar 
potential self-energy at one loop order is given by 



7^ / 77r^Tr[7oG'o((7o,q)7oGo(fco + go,k + q)] 



(27r) J (27r)3 

dqo f' d^q go(fco + go)-w2q-(k + q)-(m + Bg2)(m + B(k + q)2) 



(21) 



(27r) J (27r)3 [q^ + v^q^ + {m + Bq2)2] [(fco + qo)^ + v^k + q)2 + (m + B(k + q)2)2] 



(22) 



After some simple integrations we can show 11(0, 0) ~ and to the quadratic order in momentum and frequency we obtain 



2a;k2 1 + |(?7i2 + s2A4)i,-2A-2 



n(fco,k) 



mBv 



[l + (m + SA2)2t,-2A- 



212 



(23) 



The scalar potential self-energy causes a field renormalization of the scalar field Lp, and which in turn leads to the charge 
renormalization. 

Renormalization constants for v, m, B and a. : At this point we introduce dimensionless parameters m mv~^hr^, 
B — > Bh.v^^, and (AaA)/(27r2i;2) After collecting all the self-energies, and the Coulomb vertex corrections, and 

performing the anisotropic rescaling of the space-time coordinates as, xq — > x^e^'^, x — > xe', the parts of action without the 
disorder vertex become 



S 



(Av + Am + Ao5 + 3Aso)/ 
1 + (m + B)2 



70 (e ^'9o +ig^a) +ve ' 1 + 



2al 



STr^/l + (m + B)2 



(Ay+AM-Ao5-3Aso))^ , 
+vA I m ^ ■ — [jn + B) 



od 



-vA-M B 



1 + (to + 5)2 
al 



TT^l + {m + BY 



2,71^1 + [m + Bf 



(to + B)) X e-^'al ^tA, + 



1 



-21 



{m + B)j 

2al 1 + |(to2 + B'^)+mB^ 
[l + (m + B)2]5 



(24) 



Now we introduce the renormalization constants for the fermionic field the scalar field Lp, the charge g, the Fermi velocity 



w, the mass parameter to and the higher derivative parameter B according to ^/^ — > Z,^^ ' ip, ip ^ Z^p^^'^Lp, g ^ Zg ^'"g, 
V — > Z^^v, a — > Zg^Z^a = Z^^a, m Z^^m and B Z^^B, to recast the action in the original form, and find 



-1/2, 



-1/2, 



Zip 

z„ 



„3l 



=.(2-1)' 



(Ay + Am + Ao5 + 3Aso)l 



l + {m + BY 

2^ 



2al 1 + |(m2 +^2) + toB' 



[1 + (to + B)2]5 
2aa + |(TO2+B2)+TOB 

371- [1 + (to + B)2]5 

2a^ (Ay + Am -f Ao5 + 3Aso)' 

Sttv/H- (to + B)2 

2a/ f 



1 + (to + B)2 

al (Ay + Am - Aos - 3Aso)^ 



Zh = 



1 - 



37r^l + (to + B)2 [7r^l + (TO + B)2 
2a/ a/ 



1 + [m + BY 



B 

TO 



37r^l + (m + B)2 7rv/l + (TO + B)2 



1 + ^ 

TO 



(25) 
(26) 

(27) 

(28) 

(29) 

(30) 



Using the renormalization constants given above, we obtain the RG flow equations for v, to, B and a. Next we consider the 
renormalization of the disorder vertices due to the interplay of disorder and Coulomb interaction. The relevant diagrams are 
shown in Fig. [3jd), Fig. [3je), Fig. Of), Fig. lUg), Fig. Oh). For the renormahzation of the disorder vertices we can set all the 
external frequencies and momenta to zero in these Feynman diagrams. 
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Renormalization of Av : The renormalization of Ay from graphs Fig. Eld), Fig. Ele), arise due to the interplay of the 
pairs (l4x4,7o7j) and (70,14x4)- After accounting for the symmetry factor of two for each graph, we find the total contribution 



5A 



l(d) + l(e) _ 



4AmAsoA? 



(27r2w2)(l + (m + BA2)2i,-2A-2)2 (27r2i,2)(l + (m + BA2)2t,-2A-2)2 



(to + BA'^ 



(31) 



For the diagrams corresponding to Fig.[3jf), we take Fa = 70, and Ti, = 70, l4x4, 7o75, 7o7i, and after accounting for a 
symmetry factor of two, we find 



<5A|i^' = 2[Al + Ay Am + AyAos + SAyA^o] 



Al 



(27r2w2)(l + (m + BA2)2„-2A-2) 



(32) 



The graphs in Fig.[3jg), and Fig.[3|h) represent the renormalization of disorder couplings due to Coulomb interaction. Notice 
that Fig.[3jh) that represents a screening effect due to the presence of internal fermion loop, occurs only for random chemical 
potential. In Fig.[3jg) we set Ta = 70, and find 



SA 



Ha) 



2g^Av 



dqo f (fq 



% ~ ^ q ^ ("^ + -Sq ) 



2\2 



(27r) J (27r)3 {q^ + v^(f + {mvA + BwA-iq2)2)2q2 
After accounting for the symmetry factor of two. Fig. [3jh) leads to 







<5Af) 



4a[l + |(TO2 + B2A4)„-2A-2 



mBv 



[l + (m + BA2)2i,-2A 



-2 A-212 



(33) 



(34) 



After collecting all the perturbative corrections to Ay, we switch to dimensionless couplings, and introduce the renormaUzation 
constant We find 



^ 2(Ay + Ao5 +3Ago)? 
1 + (to + B)2 



2AmI 



(1 + (to + 5)2)2 



1 + 3(to + B)2 



2A 



so 



4al 



Ay J Sttv/I + (to + 5)2 



3^ [1 + (to, + 5)2] 2 



(35) 



Renormalization of Am : The renormalization of Am from graphs Fig. [3jd), Fig. He), arise due to the interplay of pairs 
(70j 7o)> (14x4,^4x4), (7075.7075). (7o7j.7o7i)> and (7o,7o7j)- After accounting for the symmetry factor of two for each graph, 
we find the total contribution 



^A 



l(rf) + l(e) 
M 



4(A2, + A^ + A^g + 3A|o)A; / (to + 5A^ 



(27r2w2)(l + (to + 5A2)2t;-2A-2)2 



:2\2 



y2J^2 



AAvAso^l 



(27r2i;2)(l + (to + 5A2)2t,-2A-2)2 



(36) 



For the diagrams corresponding to Fig.[3jf), we take = 14x4. and = 70, l4x4, 7o75i 7o7i' and after accounting for a 
symmetry factor of two, we find 



5A 



Hf) 

M 



2{-Ali - Am Ay + AmAos + 3Aa/Aso) 



Al 



(27r2i;2)(l + (to + 5A2)2„-2A-2)2 

In Fig.|3lg) we set Fq = l4x4, and accounting for the symmetry factor of two we find 



jm + BA ^)^ 

y2/^2 



1 



SA 



1(3) 
M 



2aAMl 



7r(l + (TO + 5)2) 



2\2 



(37) 



(38) 



After collecting all the perturbative corrections to Am, we switch to dimensionless coupUngs, and introduce the renormahzation 
constant Zam- We find 



1 - 2Ami '-'^-\^L - 2(A. - Ao5 - 3Aso)l ^ " ^"^ 



(1 + (to + 5)2)^ 

^ 4AyAgo^ 

7r(l + (j7i + 5)2)3 (1 + (to + 5)2)^ 



2al 



-A(Ai + A 



05 



3A^o)^ 



(1 + (to + 5)2)' Sny^l + (to, + 5)2 
(to + 5)2 



(1 + (to + 5)2)' 



(39) 
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Renormalization of Aqs : The renormalization of A05 from graphs Fig.|3ld), Fig.[3e), arise due to the interplay of pairs 
(7o7m lolj with i 7^ j) and (l4x 4,7075)- After accounting for the symmetry factor of two for each graph, we find the total 
contribution 



l(rf) + l(e) 
05 



4A|oAZ 



4A05AA/A/ 



(27r2w2)(l + (m + BA2)2„-2A-2)2 (27r2i;2)(l + (TO + BA2)2t,-2A-2)2 



(40) 



For the diagrams corresponding to Fig. Of), we take Fa ~ 7o75> ™d Fh = 70, l4x4, 7o75, 7o7i, and after accounting for a 
symmetry factor of two, we find 



^ . _ 2 [A^s + ApsAy - AqsAm] A^ / _ (m + BA^) 



(27r2i;2)(i + (m + BA2)2t,-2 A-2)2 



i;2A2 



6A05AS0A/ 



(27r2w2)(l + (to + BA2)2i;-2A-2)2 



(41) 



In Fig.|3lg) we set Ta = 7075, and accounting for the symmetry factor of two we find 



5A, 



1(9) 

05 



2al 



7r(l + (m + B)2)^ 



(m + BA'^ 



(42) 



After collecting all the perturbative corrections to A05, we switch to dimensionless couplings, and introduce the renormalization 
constant ^Aos ■ We find 



l + 2(Ao5 + Ay)/^^ ^ ^:^-2Aa/? ^ ' 



Aal 



(l + (m + 5)2)2 
2a?(m + B)2 



37rVl + (TO + B)2 7r(l + (m + -B)2)2 A05 (1 + (m + B)2)" 



(43) 



Renormalization of Aso '• The renormalization of A50 from graphs Fig.^d), Fig.^e), arise due to the interplay of pairs 
(7o75j7o7i)j (70j 14x4), (7o7ii l4x4)- These diagrams possess linear UV divergence, and after accounting for the symmetry 
factor of two for each graph, we find the total contribution 



°^SO 



4(2AsoAo5 + AyAA/)AZ 



4AsoAmA/ 



3(27r2w2)(l + (m + BA2)2i;-2A-2)2 (27r2w2)(l + (m + SA2)2w-2A-2)2 



(to + BA^' 



(44) 



For the diagrams corresponding to Fig.|3jf), we take Ta = Joji, and Ff, = 70, l4x4, 7o75, 7o7j. and after accounting for a 
symmetry factor of two, we find 



SA 



so 



3(27r2w2)(l + + 5A2)2v-2A-2)2 



2 [A|o - AspAy + AsqAm] Al / _ 3{m + BAy 



1)2 A2 



2Ao5AsoAZ 



3(27r2t;2)(l + (m + SA2)2„-2 A-2)2 



(45) 



In Fig.|3lg) we set Fq = 7o7i, and accounting for the symmetry factor of two we find 



2aAsol 



3n{l + {m + Byy- 



1 + 3 



(to + SA^)^ 

^^2^2 



(46) 



After collecting all the perturbative corrections to Ago, we switch to dimensionless couplings, and introduce the renormalization 
constant Zaso- We find 



Zaso = e 



-I 



, 21,. . , 1-3(to + B)2 2Am; 1-9(to + B)2 

-L - -7(^50 - ^v) 



2Ao5Z 



3' '■'(1 + (to + B)2)2 3 (1 + (to + B)2)2 l + (m + 5)2 



Aal 



AAvAmI 



2a(l + 3(to + B)2)/ 
377 (1 + (to + B)2)i 37rVl + (to + -B)2 ^ SAso (1 + (to + B)^' 



(47) 
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RG flow equations: Using the renormalization constants found for the dimensionless coupling constants we obtain the 
following RG flow equations 



dv 
dl 

dm 

~di 

dB 

H 

da 
H ' 



2a 



(A;. 



Am + A, 



05 



3Aso) 



— J— = m + (m + 3B) 



B + {m - B) 



1 + (to + 5)2 

(Ay + Aa/ - Ao, 



37r^l + (TO + B)^ 



- (m + B) 



3Aso) 



37r^l + (m + B)2 

(Ay +Am + Ao5+3Aso) 



2a 



1 + (to + Bf 



2al + |(TO2+S2)+TOB 



= Av 



- 1 



1 + (to + B)2 

2(Ay 



37r^l + (TO + S)2 



37r 



^05 



3A 



so) 



1 + (771 + Bf 



2A 



A/ 



1 + 3(to + B)2 



[1 + (to + 5)2] 5 

4q; 



(48) 
(49) 
(50) 

(51) 



{l + {m + Bfy 37r^l + [m + Bf 



4a 1 + |(to2 



B' 



mB 



[i + (to + S)2]^ 

dAM 



4Aa/Aso 

(1+ (to + 5)2)2 



(52) 



dl 



Am 
2a 



1-3(to + B)2 ^ ^ , 1-(to + B)2 

1-2Aa/-^ ^ ^ - 2(Ay - Ao5 - 3A50) ^ ' 



4a 



7r(l + (TO + B)2)2 

dAo5 



+ 



(1 + (to + 5)2)2 

4AyAso ^ 
(1 + (to + 5)2)2 



-4(A2, + A25+3A 



so) 



(1 + (to + 5)2)2 3;rv/l + (to + 5)2 
(to + 5)2 

(1 + (to + 5)2)2 



(53) 



dl 



Ao5 

4a 



6A 



so 



(1 + (to + 5)2)2 

2a(TO + 5)2 



(1 + (to + 5)2)2 (l + (m + 5)2)' 



37rVl + (TO + 5)2 7r(l + (TO + 5)2) 
dAso 



4A|o 



(1 + (to + 5)2)2 



(54) 



dl 



^so 



I -If A -A ) 1-3(to + 5)2 2Aa/ 1-9(to + 5)2 

3^ ^ (1 + (771 + 5)2)2 3 + + 5)2)2 

2a(l + 3(TO + 5)2) 



2An 



1 + (777 + Bf 



4a 



TT 



(1 + (to + 5)2)5 37r^l + (TO + 5)2 



4AyA 



3(1 + (to + 5)2)2 



(55) 



The RG equations presented in the main text are obtained by setting Am = A05 = Ago = 0. Notice that in the presence of a 
finite mass or 5, a random mass is always generated from the other disorders at quadratic order 

Fixed point analysis : The RG flow equations for generic time reversal symmetric disorder have following five fixed points 

(56) 
(57) 

(58) 

(59) 

(60) 

The fixed points FP2, FP3, FP4, FP5 describe possible universality classes of the SM-massless CDM phase transitions. 

FPl is the noninteracting, clean fixed point. Upto a critical strength of disorder this fixed point is stable, and the fermion 
mass is the only relevant perturbation. In its vicinity Aa ~ Aaoe^'- To see how the Tl-Bl phase boundary is shifted by various 
disorders, we first consider the noninteracting problem by setting a = 0. For 5, we have 5 = Bqc^K The flow equation for 
mass TO can be approximated as 



FPl 


A* 
^y 




= Ao5 = A}o =a* =m^ B 


= 0, 


z = l 






FP2 


A* 
^y 


+ AS5 


= \, ^li = A^o = a* = TO = 


5 = 


0, z = 


3 
2 




FP3 


A* 

^05 


A* 


= i At, = A*so = a* = TO = 


5 = 


0, z = 


3 

2 + 


2AX, 


FP4 


A* 


9 

~ 10' 


6 

Ago = -, Ay = AIj = a* = 


TO = 


5 = 0, 


z = 


31 
10 


FP5 


A* 

^05 


4a* 
~ "3^ 


= 1, At, = a;, = A*so = m 


= 5 


= 0, z = 


3 

" 2 



dm 

IT 



TO - 5o(A 



vo 



A 



05,0 



3Aso,o)e" 



(61) 
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which has the solution 



B 



m 



(A. 



An 



3A50) 



B 



mo 



(A 



vo 



A, 



05.0 



3Aso,o) 



(62) 



Am > Ao5 + 3 Aso 1 TI has a larger regime of stability. 

v+Am-Ao5-3Aso) 



The phase boundary shifts to m = -j(A\/ + Aj\/ — A05 — 3Aso)- If A 

In the presence of Coulomb interaction we find a further shift of the phase boundary to m = ^(Ay + Aj\/ — A05 — 3A50) ~ "§7- 
In the vicinity of this fixed point a is marginally irrelevant, and decreases logarithmically. 

The line of fixed points FP2, describes the phase boundary between the SM and massless CDM phases for the noninteracting 
problem with chiral symmetric disorder. The associated RG flow and phase diagram in Ay — A05 plane, for m = B ^ a ~ 
Am = Aso = are respectively shown in Fig. |4(a)| and Fig. |4(b)| Now we linearize the recursion relations in the vicinity of 



1.0 



0.8 



Him 




FIG. 4: (a) The RG flow, and (b) the phase diagram in Ay 
corresponds to the SM-CDM phase boundary. 



Ao5 plane for m — B 



0. The blue line Av + Ao 



this line of fixed points, and obtain 



d6a 

~dr 

dSAM 

dl 
dSAso 
dl 
dSAv 

dl 
dSAo5 

dl 



= -Sa ^ Sa = fJane'/^ 
2 

= -4A;,(5Am + 4A;,5Aso 
4 4 

^ --A*y6AsO + ^A*ydAM 

= 2A^{SAv + SAo5 + SAm ' 
= 2A*^{SAv + SAo5-5Am 



3A,o-^) 



3A 



so 



(63) 
(64) 
(65) 
(66) 
(67) 



From Eq.|64]and Eq.|65] we find 5 Am + 'iSAso ~ + 3(5A|gQ. Therefore, chiral symmetry breaking perturbations in the 
diffusive phase will remain finite. After adding both sides of Eq.l66land Eq.|62]we find 

5Av + 6A05 + 2{5Am + 3<5A5o)(At, - A^^) - ^(2A^ + A^g) - [SA^. + SA^^^^ + 2{6A\i + i5A%o){^*v ~ ^k) 



(2A^+A*5)]e' 



(68) 



is the most relevant variable with eigenvalue one. Therefore the mean free path diverges with exponent = 1. In the vicinity 
of this fixed point, the relevant variable also describes the phase boundary in the multidimensional coupling constant space, and 
z changes continuously along the phase boundary. First we note that interaction shifts the phase boundary to higher values of 



12 




Av and A 



05 



For Ay 



> Aq5 the mass and the spin orbit disorders shift the phase boundary to smaller values of Ay and A 



For Ay < Aq5 the mass and the spin orbit disorders shift the phase boundary to higher values of Ay and A05. 



05- 

Now consider 



^ -B 



dB 

m 

dm 
IF 



2a: 



Bi^2Al 



(69) 
(70) 



Notice that apart from a redefinition of the effective mass, the scaUng dimension of m has changed into (3/2 — 2Ay). This has 
important role in governing the insulator-CDM phase boundaries. For Agg > Ay, the scaling dimension of m is bigger than 
unity, and this increases the stability of the insulating phases with respect to the CDM. The RG flow in the m — Ay plane for 
B = a = Am = Aos ~ Aso = is shown in Fig.|5] 

The line of fixed points FP3 describes the phase boundary between the semimetal, and massless CDM phase for the 
noninteracting problem in the Aj\/ — Aqs plane. The associated RG flow and phase diagram in Am — Aqs plane, for 
m ^ B ~ a — Ay = Aso = are respectively shown in Fig. |6(a)| and Fig. |6(b)| The mass disorder shifts the phase 
boundary to higher values of A05, and the diffusive phase does not have chiral symmetry. Now linearizing about this line of 
fixed points we find 



dSa 

~dr 

dSAv 

dl 
dSAso 

dl 
dSAM 

dl 
d5Ao5 

dl 



4A*/Ay + AAl^SAso 
^A*/Aso + ^KiSAv 

Set 

2AIj{-SAm + SA05 - SAv + 3Aso + ^) 

OTT 



2A*^{5Ao5-SAm + SAv-3A 



so 



25a, 



(71) 
(72) 
(73) 
(74) 
(75) 



FromEq.|72]andEq.|7l we find SAy -iSAgo = (5A^ - 3(5A° . After adding both sides of Eq.|74]from Eq.|75] we find 



(5Ao5 - 6 Am + 2{5Av - 3(SAso)(Ao5 - Am) 



A6a (1 + 3A^^) 
3^ (-1+4A^,) 



5Al, - 5Al, 



2(<5A°y-3<5A"5o)(Ao5-A^ 



Ml 



A5c 



(1 + 3AXf) 



3^ (-1 + 4A*,) 



(76) 
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FIG. 6; (a) The RG flow, and (b)the phase diagram in A m 
corresponds to the SM-massless CDM phase boundary. 



Ao5 plane for m = B 



0. The blue line Aos 



is the most relevant variable if A^^ < 3/16. This variable defines the phase boundary, and critical properties are non-universal. 
For A^^ < 3/16, v = 1. For A^j > 3/16, (5Ay and provide stronger perturbation with eigenvalue 16A^^/3 > 1. The 

linearized equations for finite B and m are 

dB 

-JT = -B (77) 
al 

dm 3 1 

^ = 2"^ +2^ ^''^ 

Notice that the scaling dimension of m is 3/2, and this increases the stability of the insulating phases in the vicinity of FP3. 

The fixed point FP4 is the critical point in the class of spin-orbit disorder The associated RG flow and phase diagram in 
Aso ~ ^05 plane, for m = B = a ~ Ay = Ajv/ = are respectively shown in Fig. |7(a)| and Fig. |7(b)| Linearizing about this 
critical point we find i/ = 1. The irrelevant variable provides the phase boundary. This critical point is highly unstable against 
random potential, and mass disorder, and Coulomb interaction. The fixed point FP5 is the only finite interaction, dirty critical 
point. The associated RG flow and phase diagram in a — A05 plane, for m = B = Av = Am = Aso = are respectively 
shown in Fig. |8(a)| and Fig. |8(b)| The Coulomb interaction shifts the phase boundary to higher values of A05. Linearizing the 
flow equations we find SAy — 6Aye^\ SAm — (JA^^e^'/^, and SAso = (^^so^ From the equations for A05 and a we 
find the relevant and irrelevant combinations have eigenvalues (1 ± \/5) /2. The irrelevant variable provides the phase boundary. 



n. QUANTUM CRITICAL SCALING PROPERTIES OF 3+1 -DIMENSIONAL MASSLESS DIRAC FERMION 

In this section we consider the quantum critical scaling properties of massless Dirac fermions when disorder is irrelevant, 
and for simplicity only consider the chemical potential disorder. At a finite temperature or a chemical potential the flow toward 
infrared limit will be truncated by the largest energy scale. Thus £ = minj^T-, i^^} acts as an infra-red cutoff, where £t ~ v/T, 
and £^ ~ v/i^l are respectively the thermal de Broglie wavelength and the inter-particle separation. For £t < £^, we observe 
critical properties of the Dirac fermions, and the conventional Fermi liquid behavior is observed for £t > £fi- For this reason we 
focus on the regime £t < £fi- With decreasing temperature, the disorder strength decreases according to Ay(T) ^ AvqT/Tq, 
where the ultra-violet scale Tq = hvA/ks ^ W^K. In this regime there is an initial enhancement of a (down to a scale 
T( Ayo) < To), followed by the logarithmic decrease of a (below T( Ayo)). This leads to non-monotonic temperature dependent 
corrections to the scaling properties of non-interacting Dirac fermions. 

Specific heat and compressibility: For the noninteracting clean problem, the density of states at the Fermi point vanishes 
quadratically, p{E) oc and the free energy density has the power law dependence / ~ £-{d+z}^ j^ie free energy density is 
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FIG. 7: (a) The RG flow, and (b)the phase diagram in Ago — Aos plane for m — B — a = Av ~ Am = 0. The blue dotted line, and the 
phase boundary in (b) respectively correspond to the relevant and irrelevant variables at FP4. 
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FIG. 8: (a) The RG flow, and (b)the phase diagram in q — Aos plane for m = B = Av ~ Am = Ago = 0. The blue dotted line, and the 
phase boundary in (b) respectively correspond to the relevant and irrelevant variables at FP5. 



given by 



Li4(-e^/'^) +Li4(-e-''/'^) 



(79) 



In the limit /x/T < 1, 



f{T,fi)^ 5- 



180 6r2 127r2T4 



(80) 
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The specific heat follows the scaUng relation C ^ T'^^^, and is given by 

C 



-T 



QT2 



2^3 



V 



15 



J£_ 

3T2 



The compressibility also follows expected scaling relation k 



rpd/z 1 ^jjj given by. 



In the Hmit fi/T <^ IWe find k 



,-3 



,2 



Li2(- 



)+Li2(- 



(81) 



(82) 



+ f-i^/n^). For /.j = 0, the ratio C/{kT) = Tir^kj^/b is a universal number 



Due to interplay of disorder and interaction, there is an initial enhancement of C, and k followed by the logarithmic suppression 
by a factor (1 + 4ao/(37r) \og{To/T))^^^^ . The non-monotonic behavior becomes pronounced if Ayo is close to the critical 
strength. 

Diamagnetism: We first consider the clean, noninteracting problem. The magnetic field B introduces another length £b = 
(eB)-i/2. AtT = /i = 0, a naive application of the scaling formula gives / ^ (eS)^, and a constant diamagnetic susceptibility 
X- In the presence of gauge field naive scaling hypothesis becomes inapplicable. The interaction of electrons with an external 
gauge field requires the use of RG scheme with a proper regularization procedure to address the ultraviolet divergence of the 
fermionic polarization bubble. Since d+ z ~ 4, the problem is at upper critical dimension and one should anticipate logarithmic 
corrections 111]. A proper analysis leads to the renormalization of charge and the gauge field, and a (eB)^ log(a/£) contribution 
to the free energy density ||2l|3|]. Again £ has to be chosen to be the smallest of the scales chosen from (ei?)^^/^, v/T, v/ii, 
and v/m, and accounts for the large value of x- Such logarithmic enhancement has been argued for the large diamagnetic 
susceptibility of Bi, and the narrow gap semiconductors such as Pbi_i:Sn2;Te, Bii_a;Sbx. 

The Landau level spectrum of massless Dirac fermion is En.k, = ±\/v^k^+E^n; where Eb ~ \f2vllB is, and n is the 
Landau level index. The degeneracy of n = level per unit area is eB/2'K and that for levels with n > is eBjix. The 
diamagnetic susceptibility arises from n 7^ levels, and in the following we consider the part of the free energy density E that 
arises from n 7^ levels. This is given by 

r(-l + |) / 



E = lim 



I V2 



"5 7724 



The divergent part of E is given by 



dx 



l/2-t/2 



n=l 



lim 

e-i-O 



V2j r(^ + f) 



(83) 



E, 



and it is absorbed into the vacuum energy density 5^/2, which leads to the field and the charge renormalizations 



Bi 



B^ [I 



(84) 



(85) 



(86) 



In the above equations 1/e ~ log(A£B/V2), and we have restored the explicit dependence on the speed of light c. The finite 
part of the energy density is 



E 



finite 



C(-l)log 



+ C'(-l) + C(-l)W2)-^(-l/2)) 



(87) 



where ip is the digamma function, and ( is the Riemann zeta function. From the finite part of the energy density we find the 
diamagnetic susceptibihty 



X ' 



(log f ^ 
247r2 B 



1.74 ) , ge < V2£t 



(88) 



where Bq ~ hh? /{2e) ~ lO^T. At a finite temperature, a similar calculation can be performed for the free energy density and 
in the high temperature limit we find 



24^2 I y 



1.74 1 , Ib > V2in 



(89) 



The departure from the noninteracting formula follows from the renormalization of e?v oc av^. When a decreases logarithmi- 
cally, we can set z ~ 1, and find e^v = constant. Therefore Coulomb interaction does not modify the scaling behavior of x- 
The disorder only causes a small suppression of noninteracting value of x- 
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Dynamic conductivity: Consider again the noninteracting, clean limit. Since conductivity a ^ t and 2; = 1, we can 

write G{bj, T) = e^T / {2ttv)<^{ijj /T, fi/T), where the scaling function 



4>(x,y) = l/36[(87r2 + 2V)^(a;) +3x{tanh(2;/4 + y/2) +tanh(.T/4-y/2)}] 



(90) 



Again we will focus on ^ ~ 0, or (a;, T)^ ^ limit. The inelastic scattering rate ^ a^T, is larger than the elastic scattering 
rate due to disorder, and the conductivity will be mainly governed by the interaction effects. If w ^ a^T, the leading order 
answer for the dynamic conductivity follows from the noninteracting formula, 



12ttv 



tanh(w/4T) 



(91) 



As <j{uj, T) cx a, there is an initial enhancement, followed by a logarithmic suppression by the factor a [l 



4q, 
Sir 



aloe 



In the opposite limit w <C a^T the collision processes in the particle-hole plasma due to Coulomb interaction governs the 
conductivity iflBt]. To obtain concrete answer we have performed a calculation using the quantum Boltzman equation 



(9t + eE- Vp)/a(p,t) =-C[/„](p,t) 



(92) 



within leading log approximation (|4l|51]. In the above equation E is the external electric field and /a(p,i) is the fermion 
distribution function, and a is the collective label for particles and holes, and also the chiralities. The definition of the collision 
operator C[/a](p,t) involves the square of the amplitudes of the two particle scattering processes (particle-particle, hole-hole 
and particle-hole), and a combination of Fermi functions. In order to solve the Boltzman equation we introduce the ansatz 



/a(P,^) 
/a°(P) 



P 



2^5{lo) + (1 - /,"(p))e,E(c.) • ^Xa(p,c. 



(93) 



where is equilibrium distribution function, and convert the linearized Boltzman equation into a varitaional problem for 
Xa(p, uj) which are functions of dimensionless variables |p|/r andw/T. In the particle-hole symmetric case /i = 0, x+(Pi '-^) = 
— X_ (p, Lu), and particles and holes equally contribute to the transport. At the end we extremize the functional Q[x], given by 



Q[x\ 

2^2 



dpf{p)(l-f{p)) 



1447^3 



2T' 



(94) 



Now choosing a single parameter ansatz x{Pi = {p/T)"g{u!/T), we find the extremum occurs for n ~ 0.896, and 



a{LO,T) 



30.46r 



iuj 26.67 



(95) 



alog{l/a) 

One popular choice n — 1 only accounts for particle-hole scattering. The proximity of our n ^ 0.896 to 1, suggests that the 
contribution from particle-particle collisions is small. Only for fi > T, the like-particle collisions will be dominant and lead to 
conventional Fermi liquid result ^ a'^T'^/ fi, and ctq ^ M'^/(Q!T^)- From the expression for dynamic conductivity we can 
see the existence of a Drude peak, and T-linear dc conductivity. The renormalization of a will now cause an initial suppression, 
followed by a logarithmic enhancement of the dc conductivity. 
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